Evolution of Edge states and Critical Phenomena in the Rashba Superconductor with Magnetization 
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We study Andreev bound states (ABS) and resulting charge transport of Rashba superconductor (RSC) where 
two-dimensional semiconductor (2DSM) heterostructures is sandwiched by spin-singlet s-wave superconductor 
and ferromagnet insulator. ABS becomes a chiral Majorana edge mode in topological phase (TP). We clarify 
that two types of quantum criticality about the topological change of ABS near a quantum critical point (QCP), 
whether ABS exists at QCP or not. In the former type, ABS has a energy gap and does not cross at zero energy 
in non-topological phase (NTP). These complex properties can be detected by tunneling conductance between 
normal metal / RSC junctions. 
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Topological quantum phenomena and relevant quantum 
critic alityhave been an important concept in condensed matter 
physics Recently, stimulated by the issue of Majorana 

fermion in condensed matter physics yj-lg], topological quan- 
tum behavior of superconductivity becomes a hot topic 
1211 . One of the most crucial point is the property of the non- 
trivial edge modes in topological phase where edge modes are 
protected by the bulk energy gap. 

The edge state of superconductor has been known from the 
study of Andreev bound state (ABS) in unconventional su- 



perconductors II 13141 31 . In high Tq cuprate, dispersionless 



zero energy ABS ubiquitously appears IU41I15H due to the sign 
change of the pair potential on the Fermi surface. The zero en- 
ergy state manifests itself as a zero bias conductance peak in 
tunneling spectroscopy II15[ Il6ll . Subsequently, the presence 



of ABS with linear dispersion has been clarified in chiral p- 
wave superconductor II 1 711 realized in S^RuCU, where time 



reversal symmetry is broken 11811 . On the other hand, in the 
presence of spin-orbit(SO) coupling with time reversal sym- 
metry, it has been revealed that spin-singlet s-wave pairing 
and spin-triplet p-wave one can mix each other due to the bro- 
ken inversion symmetry |T^-21|. ABS appears as a helical 
edge mode appears for A p > A s where we denote s-wave 
and p-wave pair pote ntials a s A s and A p , respectively, with 



A s > Oand A„ > 12 U 1221 



The critical behavior of ABS has been discussed in spin- 
triplet chiral p-wave pairing yfl. By changing the chemical 
potential pi of spin-triplet chiral p-wave superconductor from 
positive to negative, ABS as a chiral Majorana mode disap- 
pears. The corresponding quantum critical point is /t = 0. 
Although, such a quantum phase transition can be possible in 
v = 5/2 fractional quantum Hall system and cold atom 
JHUl, it is significantly difficult to obtain superconducting 
state for negative p in electronic superconductors. 

In all of above works, ABS is generated from unconven- 
tional pairing with non-zero angular momentum. On the other 
hand, in the presence of strong SO coupling with broken time 
reversal symmetry, chiral Majorana modes can be generated 
from spin-singlet s-wave pairing J25, 26]. Fu and Kane have 
revealed the presence of chiral Majorana mode at the bound- 
ary between ferromagnet and superconductor generated on the 



surface of topological insulator (TI). After that manipulating 
Majorana mode in TI [26J| and in semiconductor hetero struc- 
tures based on conventional spin-singlet s-wave superconduc- 
tor have been proposed in several contexts B27H29I1 . Sau et 
al. has proposed a unique Rashba superconductor where two- 
dimensional electron gas (2DEG) is sandwiched by conven- 
tional spin-singlet s-wave superconductor and ferromagnetic 
insulator 12811 . These systems are really promising for future 
application of quantum qubit since host superconductor is ro- 
bust against impurity scattering. 

Although there have been several theoretical studies about 



the present RSC O0M32I1 . the feature of the Andreev bound 



state (ABS) and its relevance to the topological quantum 
phase transition has not been revealed at all. It is known that 
ABS emerges as a chiral Majorana edge mode in TP, however, 
the evolution of ABS in the non-topological phase (NTP) and 
its connection to quantum phase transition have not been clar- 
ified yet. To reveal these problems is indispensable to under- 
stand the tunneling spectroscopy of normal metal /RSC junc- 
tion system and future applications of quantum device. 

In this Letter, we study energy dispersions of ABS in RSC 
composed of 2DEG sandwiched by spin-singlet s-wave su- 
perconductor and ferromagnetic insulator. It is clarified that 
there are two types of quantum criticality for ABS, i.e., quan- 
tum phase transition with or without ABS corresponding to 
type I and type II, respectively. In type I, ABS can exist 
even at critical point where bulk energy gap closes and in 
the NTP. Nonzero ABS generated in the NTP does not cross 
at zero energy. These features are completely different from 
those in type II where edge states become absent both at the 
critical point and in the NTP. The conventional criticality of 
spinless s pin- triplet chiral p-wave superconductor belongs to 
type II I3I. I24I1 . The conductance between normal metal / RSC 
junction shows wide variety of line shapes reflecting on these 
novel quantum criticalities. We also show the drastic jump of 
the conductance at critical point. 

A Hamiltonian of Rashba superconductor with magnetiza- 
tion is given by the following form J^-Hit] : 

H(k) = H (k)+H R (k) + H z + H s , (1) 

where kinetic energy Hq, Rashba spin-orbit interaction 
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(a)mA 2 <|K| (b) m.A 2 > |K| 
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FIG. 1. (color online) Energy spectra of the normal (A = 0) states. 
(a)Zeeman (Rashba spin-orbit) interaction is dominant with m\ 2 < 
\V Z \. (b)Rashba spin-orbit interaction is dominant with m A 2 > \V Z \. 
The critical value of chemical potential for the transition between 
topological and non-topological superconductors is given by ±/x c = 
±s/V£ — A 2 , (see discussion below eq. 




FIG. 2. (color online) Normal metal (N) / Rashba superconductor 
(RSC) junction. Andreev bound state as edge state can exist denoted 
by the (red) arrow. 



(RSOI) Hr, Zeeman interaction Hz by exchange field from 
FM insulator, and spin-singlet s-wave pair potential H§ in- 
duced by proximity effect are Ho(k) = £feSoT z , H^(k) = 
X(s x T k y - s y T z k x ), H z = V z s z t z ,H s = -As y T y , where s 
and r are Pauli matrices, sq and tq are 2 x 2 unit matrices, de- 
scribing electron spin and particle-hole degrees of freedom, 
respectively. We take the explicit form of kinetic energy as 
£fc = k 2 /2m — fi with fx being chemical potential, for simplic- 
ity. The exchange energy in a 2DEG can be tuned by changing 
the material of ferromagnetic insulator, or tuning the barrier 
thickness between the ferromagnetic insulator and the 2DEG. 
In the normal states ( A = 0), there are two types of the energy 
bands as shown in Fig. Q] For Zeeman interaction dominant 
case with mX 2 < \V Z \, there are two parabolic dispersions 
(Fig- Q2 a ))- On the other hand, for RSOI dominant case with 
m\ 2 > \V Z \, the shape of the energy band is wine-bottle like 
(Fig. |TJb)). As we shall see later, the difference between these 
two types of energy bands in normal state becomes important. 

The eigenvalues of the Hamiltonian for the infinite sys- 
tem are_given by E a (k x , k y ) = s/r)k + Ck, E b (k x ,k y ) = 
- Vfjk + Ck , E c {k Xl k y ) = vV - Ck, and E d (k x ,k y ) = 
-VVk - Ck with 



2;„2 



m = + yk 



a = 2^(\ 2 k 2 + v 2 )e k + v 2 A 2 , 



(2) 



k 2 with real k x and k y for 



where k is defined by k = yk 2 

the plane wave. The corresponding eigenvectors u a (k x , k y ) 
with a = a, b, c, and d are also obtained analytically. 

Let us now consider a semi-infinite RSC in x > with flat 
surface at x = 0. The wave function in the present system is 
given by 



iP ky , E (x > 0) = Y / t i u i (q i ,k y )e i * x e ik *y. (3) 

When qi is a real number, the corresponding wave function ex- 
presses propagating wave, i.e., scattering state. On the other 
hand, when qi is a complex number, it describes an evanescent 
wave. Energy E and y-component of momentum k y are good 



quantum numbers. To obtain qi, we solve k for fixed E = 
E a (k x ,k y ) and E = E c (k x ,k y ) for E > [E = E b (k x , k y ) 
and E = Ed(k x , k y ) for E < 0]. qi is given by = k x 
by postulating the constraints dE a (qi, k y )/dqi > for scat- 
tering state, and Imqt > for evanescent state. Note here 
that, in general, k and qi become complex numbers which can 
be obtained by analytical continuation. The coefficient ti is 
determined by the confinement condition as ^^(O) = 0. 

Tunneling conductance of normal metal (N) / RSC junction 
as shown in Fig. [2] is calculated based on the standard way 
Jl5 . 331. Suppose that the normal metal has no spin-orbit 
interaction, i.e., the Hamiltonian reads Huik) = (k 2 /2m — 
Hm)sqt z , where /im = H — eo with eo being the energy of 
bottom of the energy band, which is negative, and the interface 
potential is given by H\ = HsqTqS(x). The wave function in 
N is given by 



lpk y ,E,s{x < 0) 



~T" ^ ^ Tss't'Xs't' 



-it' k i x 



Jk y y 



(4) 



where the first term denotes an incident electron with spin s, 
and Xst is the eigenvector of spin s for electron (r = +1) 
or hole (t = —1), and k cx = ^/2m(/iM + E) — k 2 and 



k\ lx = y 2m(/iM — E) — k 2 are momenta of reflected elec- 
tron and hole, respectively. On the other hand, the wave func- 
tion in RSC (x > 0) obeys the same form as in eq. Q. 
The boundary condition at the interface located on x = 
is given by the following two expressions 13411 . ip(— 0) = 
VX+0), w(+0)^(+0) - w(-0)V>(-0) = -i2#T z ^(0), where 
velocity in .T-direction is v(x) = dH/dk^k^^-id^ - Solving 
the above equations, we obtain reflection (transmission) coef- 
ficient r (t). Charge conductance G normalized by its value 
Gn in the normal state (A = 0) with V z = 0, [i/m\ 2 = 4, 



Hm/mX 2 



2 x If) 4 , and Z 2 



mH 2 /fi 



M 



10 4 , which 



corresponds to the case of Figs. |3jh) and|3k) with A = 0, at 
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FIG. 3. (color online) Energy spectra and tunneling conductances as 
a function of bias voltage (eV/ A) of the Rashba superconductor. The 
horizontal axis denotes the normalized momentum k = fej/VmA. 
Zeeman interaction and Rashba spin-orbit interaction are fixed as 
Vz/A = 2, mA 2 /A = 0.5. The chemical potential is set as fol- 
lows. (a),(d):/VA = 0, (b),(e):/x/A = 1.7, (c),(f):/x/A = V%, 
(g),(j):/x/A = 1.8, (h),(k):/VA = 2, (i),(l):/x/A = 2.5. 



zero bias voltage (eV = 0) is given by 



G/G N = J2 / dk y T s (k y ,E) 
/j2 f ¥ d kyT s (k v ,0), 



(5) 



with T s {k y ,E) =2 - J2 s > t > t '\i'ss't'\ 2 and p M = 2mk 2 F . 
Hereafter, the parameters are fixed as Z 2 = 10, hm/A = 10 4 , 
and all the conductances G are normalized by the same value 
of G N . 

We discuss the energy spectra and the tunneling conduc- 
tances, focusing on the difference of the criticality between 
two RSCs with different chemical potential with \i > (Fig. 
0) and fj, < (Fig. 0) for | V z \> m\ 2 . 

In TP (Fig. |2a) and Fig. HJa)), ABS appears as a chiral 
Majorana edge mode, where \V Z \ > ^/ /j, 2 + A 2 is satisfied. 
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FIG. 4. (color online) Energy spectra (upper) and tunneling conduc- 
tances (lower) of the Rashba superconductor for negative chemical 
potentials. (a),(d):/i/A = -1, (b),(e):^/A = -y/3, (c),(f):/x/A = 
—2. The other parameters are the same as in Fig. [3] 



Due to the presence of this mode, the corresponding tunneling 
conductance has a zero bias peak as shown in Fig. [3jd) and 
Fig. Ed). For n > 0, near the QCP [Fig. Ob)], although 
ABS appears as a chiral Majorana mode, the corresponding G 
has a zero bias dip as shown in Fig. Oe) due to the presence 
of a parabolic dispersion of bulk energy spectra near k y = 0. 
At QCP (Fig. EJc)), it is noted that ABS remains although 
the bulk energy gap closes at k y = 0. This feature is quite 
different from (j, < 0, where ABS is absent at QCP (Fig. @] 
(b)). The resulting G has a I^-shaped zero energy dip both 
for two cases shown in Figs. [3{f) and He). For /i > 0, ABS 
still remains even in the NTP as shown in Figs. Hg),|3lh), and 
|3li). ABS has an energy gap and is absent around k y = 0. The 
tunneling conductance shows a gap structure around eV = 
[Fig. Elj)]. With the increase of fj,, i.e., away from QCP, the 
additional non-zero ABS around k. y = [Fig. Oh) and[3ji)] 
with the almost flat dispersion are generated. As a result, G 
has two peaks at the corresponding voltages inside the bulk 
energy gap (Fig. |2k) andQl)). On the other hand, for (j, < 0, 
ABS is absent in NTP as shown in Fig. He). The resulting G 
is almost zero inside the bulk energy gap (Fig. 2Jf)). Based on 
these results, we can classify two types of criticality whether 
edge state exists at QCP or not. We denote former type as type 

I and the latter one as type II in the following. 

We have also studied for | V z \ < m\ 2 . The energy spectra at 
QCP with positive fi (Fig. |5ja)) and negative fi (Fig. |5Jb)) are 
shown. In this case, irrespective of the value of /i, ABS exists 
at QCP. Therefore, the resulting criticality is always type I. 

Type I and II transitions can be distinguished experimen- 
tally by the line shape of G. In type I transition, line shape 
of G becomes almost symmetric with respect to eV = as 
shown in Figs. Of), 0c), and|3d) as compared to that in type 

II as shown in Fig. Uf). Furthermore, G at type I transition 
takes one order of magnitude larger value than that at type II, 
due to contribution from the edge states. 

It is noted that the small value of Z 2 does not qualitatively 
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FIG. 5. (color online) Energy spectra and tunneling conductances 
of the Rashba superconductor for mA 2 > \V Z \ at quantum critical 
point. (a),(c):/VA = \/T25, (b),(d):^/A = -VT25. The other 
parameters are taken as follows. mA 2 /A = 5, V 2 /A — 1.5. 



FIG. 6. (color online) Conductance as a function of chemical poten- 
tial /i/A and Zeeman interaction V z /A. Rashba spin-orbit interac- 
tion is taken as mA 2 /A = 0.5. The transition of type I (II) occurs 
at positive (negative) [i. The solid (broken) line indicates the critical 
line of type I (II) transition. 



change the results of the paper In the low transparency limit, 
the contribution from edge states becomes dominant for the 
conductance G, then the resulting line shape of G becomes 
insensitive to the parameters of the normal metal, i.e., Z 2 , 
Hm, and to. In the present case, the transmission probabil- 
ity in the normal state (A = 0) becomes sufficiently small 
with Gn/Go ~ 10 4 , where Go denotes the maximum value 
of Gjy, even for Z 2 = since the magnitudes of Fermi mo- 
menta in left normal metal (x < 0) and right RSC (x > 0) are 
much different with /Um/a* > 10 3 - 

Here, we mention the criticality of ABS in spinless chiral p- 
wave superconductor. Hamiltonian of spinless chiral p-wave 
superconductor is given by 



H p (k) 



k 2 /2m — fx 
A p fc+ 



A p fc- 



-k 2 /2m + fj, 



(6) 



It is known that QCP is located at fi = 0. ABS appears as a 
chiral Majorana mode in TP (/i > 0) while it is absent in NTP 
fx < 0, respectively yfl. ABS disappears at QCP. In the light 
of our classification, quantum criticality of spineless chiral p- 
wave superconductor belongs to the type II. 

To understand the difference of two types of criticality, we 
focus on the energy dispersions in the normal state shown in 
Fig. Q] Here we introduce the critical value of transition be- 
tween TP and NTP ±li c = ± V / V 2 - A 2 . The ABS is gener- 
ated from —hp to where the magnitude of kp is almost 
the same with that of the large Fermi surface. First, we focus 
on the case with to A 2 < \V Z \. The type I quantum phase tran- 
sition occurs at /i = li c , shown in Fig. [3] In this case, the 
large Fermi surface survives as shown in Fig. QI a )- On the 
other hand, as shown in Fig. [4] type II quantum phase tran- 
sition occurs at fi = —/j. c . In contrast to the type I, the large 
Fermi surface vanishes in the NTP as shown in Fig. [TJ a )- F° r 
m A 2 > \V Z \, the quantum criticality always belongs to type I. 
Actually, as shown in Fig. |TJb), the large Fermi surface sur- 
vives both at fx = (i c and fi = — /j, c . For type I, the number of 
Fermi surfaces is 2 in NTP and 1 in TP. On the other hand, for 



type II, the number of Fermi surface is in NTP and 1 in TP. 
Above rich behavior of quantum criticality in RSC originates 
from the simultaneous existence of the Rashba spin-orbit cou- 
pling and the Zeeman interaction. 

Finally, we show the zero-bias tunneling conductance of 
RSC as a function of /i and V z in Fig. [6] The quantum phase 
transition from NTP to TP occurs with tuning the parameter 
V z or fi. In accordance with this transition, the conductance 
increases by about three orders of magnitude, due to the con- 
tribution from zero energy ABS at k y = 0. 

In this letter, we have calculated the energy spectrum and 
the tunneling conductance of RSC and clarified its quantum 
criticality. Quantum phase transition between topological and 
non-topological superconductors has two types of criticality 
whether ABS survives or not at QCP. It is remarkable that 
ABS can remain at QCP in RSC distinctly from spinless chi- 
ral p-wave superconductor which is a prototype of topological 
superconductor. This stems from the structures of Fermi sur- 
faces which are spin-split by Rashba spin-orbit interaction in 
the normal state. This results can provide a new perspective of 
quantum criticality for topological superconductors. We have 
considered only the spin-singlet s-wave superconductor. It is 
interesting to study in the case of unconventional supercon- 
ductor where much richer quantum criticality can be expected 
|35rSl. 
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